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Abstract 

After a general overview of some features of the relaxation dynamics of the Hamiltonian Mean Field model, its equilibrium 
thermodynamic properties are used to rephrase the out-of-equilibrium regime for energies below the critical point u c = 0.75 
in terms of an effective dynamical coexistence between a clustered and a gaseous phases, whose existence could be associated 
to the large relaxation times observed when Ui < u < u c , with u\ = 0.5. Starting from the hypothesis that the parametric 
resonance is the microscopic mechanism allowing the energetic interchange between the particles during the collisional regime, 
a phenomenological Fokker-Planck equation based on a Langevin equation with a multiplicative noise is proposed in order to 
describe the collisional relaxation of this system towards its final equilibrium, which supports the following dependence of the 
collisional relaxation timescale r cr = tqN = IN/g. 

PACS numbers: 05.20.Gg; 05.20.-y 



I. INTRODUCTION 



II. GENERAL CHARACTERIZATION 



Despite of the Hamiltonian Mean Field (HMF) model 
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is a ferromagnetic toy model enough amenable for allow- 
ing an accurate numerical and analytical characteriza- 
tion, it exhibits many features observed in more realis- 
tic long-range interacting systems such as: violent relax- 
ation, persistence of metaequilibrium states, slow colli- 
sional relaxation, phase transition, anomalous diffusion, 
etc. These features support the reason why it can be 
considered as a paradigmatic toy model of the real long- 
range interacting systems [l4| . 

The incidence of long-range interactions does not al- 
low to divide this system into independent subsystems 
even in the thermodynamic limit. This fact evidences 
its intrinsic nonextensive nature, which distinguishes the 
HMF model from other ferromagnetic models with short- 
range interactions despite they share many analogies in 
most of the equilibrium thermodynamic properties 
The nonextensivity of the HMF model is particularly im- 
portant to understand its nontrivial dynamical behavior, 
which is still an open problem attracting much attention 
in the last years [lIL E3, EE El, d, d, E3, d] • 

The present paper of this series [l9[ is the first work 
devoted to study the dynamical behavior of the HMF 
model with a special emphasis on the relaxation processes 
undergone by this system during the evolution towards 
the thermodynamic equilibrium. We begin performing a 
general characterization of the most important features 
of the microscopic dynamics by taking into consideration 
the results obtained from the analysis of its equilibrium 
thermodynamics [l9j]. Such a study possibilities us to 
conjecture the equilibration mechanism leading this sys- 
tem to its final relaxation, which is here used to propose 
a phenomenological approach of the collisional evolution. 



A. Numerical computation of microscopic dynam- 
ics and the Vlasov equation 



According to the Hamiltonian 
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the motion equations of the HMF model are given by: 



1 



Li, Li — gN (m y cos B l - m x sin 0,) , (2) 



for i — 1,2, ...N, where m x and m y are the Cartesian 
components of the magnetization vector m = (m x , m y ) = 

(j2iLi m i) / N where m, = (cos l , sin 0j). 

From a dynamical viewpoint, the HMF model repre- 
sents the long-range interacting version of harmonic os- 
cillators system [lj] since the individual dynamics of a 
given rotator can be easily rephrased as the dynamics of 
the mathematical pendulum: 



+ msm(di -6) =0. 



(3) 



We have used the polar representation of the magneti- 
zation vector m = m (cos sin#) and the characteristic 
units for time and momentum: 




VTgN, 



(4) 



in terms of the moment of inertia I and the coupling 
constant g. Obviously, the characteristic time unit tq 



provides the timescale for the evolution of each rotator, 
and hence, tq is the characteristic microscopic timescale. 
The use of characteristic unit Lq allows to express the 



total energy per particle u = k 



1 — m ) and the 



kinetic energy per particle k = K/N — \ \pf^j /N, 

with pi = Li/Lg, in terms of the characteristic unit Sq = 
E /N = gN introduced in our previous paper [l9| . 

Another convenient way to rewrite the dynamics of 
the HMF model © follows from the introduction of the 
tridimensional vectors m; = (cos 0j, sin#i, 0) and = 
(0,0,pi): 



m, = o, x m, , o, = m, x m. 



(5) 



where the magnetization vector m obeys the following 
dynamics: 
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x m, 



(6) 



being Y the magnetization rate of change. 

A simple inspection of Eqs.Q, (O and © allows to 
understand that the dynamical phenomenology the HMF 
model can be explained from the conjugation of the dy- 
namical features of the mathematical pendulum and the 
character of the magnetization vector m evolution. Since 
the individual terms x mfc in the magnetization rate 
of change Y have a undefined signature, the dynamical 
evolution of magnetization vector crucially depends on 
the correlations among the individual rotators evolution. 
If present, the magnetization experiences large variations 
in a timescale comparable with the microscopic time To, 
where takes place a very effective energy interchange 
among the rotators. This dynamical regime character- 
ized by the existence of a collective motion completely 
analogue to the one operating in the astrophysical sys- 
tems, which is usually referred as violent relaxation [201 ] 
and very- well described by the Vlasov dynamics: 



d_ 

dt 



d_ 



m(fl)xm[/]|j/(^;t) = 0, (7) 



with m [/] = J m (9) f (8,p; t) dddp, which accounts the 
collisionless regime of the one-body distribution func- 
tion f (0,p;t). This fast relaxation regime finishes after 
arriving at certain stable quasi- stationary state (QSS) 
of Eq. Q , / (0, p; i) - f QSS (e,p) = F[e (p, 9)] , where 
e (p, 9) = |p 2 — m • m (6*) is the energy of an individual 
rotator. Hereafter, the dynamical evolution depends on 
collisional effects that decrease with the system size N, 
and hence, the system evolves throughout stable QSSs to- 
wards its final equilibration [H G3, 111, HI EE HI LLZl El • 
During this quasi-stationary regime, the individual ro- 
tators becomes almost non correlated, since the two-body 
correlation function probably survives at the 1/N ap- 
proximation as its microcanonical estimate, Eq.(58) of 



ref.[19(. Therefore, their short-time dynamics can be de- 
scribed as free motions of a mathematical pendulum. The 
vector Y experiences zero mean fluctuations in the mi- 
croscopic timescale tq whose amplitudes decrease with 
the system size as 1/yN, provoking in this way a very 
small fluctuating behavior of the magnetization vector m 
around certain mean value. Such fluctuations constitute 
the microscopic relaxation mechanism allowing the ener- 
getic interchange among the rotators, which reduces its 
effectiveness with the increasing of the system size N. 
This last observation clarifies the reason why the charac- 
teristic relaxation timescale of such a collisional regime 
T cr ~ T eq grows with N . 

The above picture is easily verified by means of micro- 
canonical numerical computation of dynamics For 
comparison purposes, we also carry out the numerical 
integration of the collisionless dynamics (O by imposing 
a truncation of distribution function f(9,p;t) = for 
\p\ > Pc and introducing a second-order finite differences 
scheme as follows: 



dt 



- PjAefij + x m [/] A p /y = 0, (8) 



where fij = f (8i,pj) with 9$ = ih\ and pj = j/12, being 
i = [l,2,...,Aq + 1], j = [-N 2 -1,...,N 2 + 1], h x = 
2ir/ (Ni + 1) and h 2 = p c / {N 2 + 1), with the boundary 
conditions: 



iVi+1,3 



fl,j 5 fi, 



±(N 2 + 1) 



0. 



(9) 



Besides, magnetization is rephrased as m [/] = 
hih 2 ^2ij nij/y with rxii = m (9i), and the partial numer- 
ical derivatives are given by Asfa = (/i+ij - fi,j) /hi 
and Apfij = (fi,j+i - fi,j-i) /2h 2 . Numerical integra- 
tions of the microscopic dynamics ([2]) and the scheme 
© with Nx = 101 and N 2 = 200 where performed by 
using fourth-order Runge-Kutta method with constant 
timesteps 5t m = 0.05 and St v — 0.02 respectively, which 
ensure a good enough conservation of the energy per par- 
ticle u. All microscopic simulations starts from water 
bag initial conditions (WB): \m(t = 0)| = 1 (Vi, 9i = 0) 
with a random uniform distribution for rotators momenta 
[ — P*,P*] with = v&u, an unstable initial condition 
which can be considered as microscopic configuration 
highly ordered. In order to avoid abrupt changes in the 
partial numerical derivatives Ag/y and A p /y , the WB 
initial conditions of the microscopic dynamics were ap- 
proximated by the following initial condition for the one- 
body distribution function: 



f wb (9,p) oc 



with Afl 



1 



1 



1 + exp [\ 9 - B*)] 1 + exp [X p (\p\ - p*)} ' 



X p = 50 and momentum cutoff p c = 1.5p*. 
The small parameter 9* used here ensures a precision of 
the initial magnetization as S |m (0)| ~ 10~ 3 . As in many 
studies 0, H, E3, HH, E3, E3] , the present calculations were 
performed at u — 0.69, an energy value corresponding at 
the equilibrium to the ferromagnetic state close to the 
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FIG. 1: Evolution of the nonequilibrium temperature To ob- 
tained from the microscopic dynamics (MD) with N = 5000 
and 10000, which can be compared to the Vlasov collisionless 
dynamics. Horizontal scale represents the function In (1 + t) 
in order to exhibit the whole simulation time interval by using 
a linear scale for small t and in logarithmic scale when t > 1. 





FIG. 2: First panel: Dynamical evolution of the nonequilib- 
rium temperature associated with two realizations of the wb 
initial conditions with N = 100000, where it could be appre- 
ciated a dynamical sensibility to the initial conditions. Panels 
(a-h) display different instants (indicated in the first panel) of 
the dynamical evolution of the system in the p-space, where it 
is revealed the formation and destruction of filamentary struc- 
tures. Black and red points are use to illustrate the mixing of 
microscopic dynamics. 



critical point u c = 0.75 of the continuous phase transi- 
tion. 

FIG Q] shows the evolution of the nonequilibrium tem- 
perature To = 2K/N (two times the kinetic energy per 
particle) obtained from two individual trajectories of ([2]) 
with u = 0.69 for N = 5000 and 10000 respectively. No- 
tice that the temporal dependences of To derived from 
the microscopic dynamics can be compared to the evolu- 
tion of this same quantity obtained from the numerical 
implementation of the Vlasov dynamics. 

Size effects are hardly appreciated at the first stages 
of the microscopic dynamics where takes place a very 
fast relaxation followed by a short transcend regime (the 
plateau in To), which finishing with another relaxation 
regime where take place the development of resonances 
(peaks) with the incidence of finite size effects. The fact 
that the immediate response of the system can be ac- 
counted for by the Vlasov dynamics demonstrates the 
existence of a collisionless evolution during this initial 
relaxation regime, whose relaxation time r vr is compa- 
rable to the microscopic timescale to, t vt ~ To. FIG[T] 
also shows that the size effects and fluctuations are very 
important in the subsequent evolution after the violent 
relaxation. The large transcend regime evidences the oc- 
currence of quasi-stationary evolution where the nonequi- 
librium temperature To remains a long time fluctuating 
around a nonequilibrium value (compare it with the equi- 
librium temperature T eq — 0.476 represented as a hori- 
zontal line). 

The evolution of a given realization of the microscopic 
dynamics in the /Lt-space (0,p) is shown in FIGJUfor N — 
100000. It illustrates how the highly ordered initial state 



is progressively destroyed by the incidence of chaoticity 
and mixing of the collisionless regime. Such mechanisms 
manifest in the development of filamentary structures, 
which disappear at a coarsed grained scale leading thus 
to the establishment of the stable QSS. This process is 
also displayed in FIG|3l but this time, the evolution of the 
one-body distribution function obtained from the finite 
differences scheme (J5J) where it could be also noticed the 
formation and destruction of filamentary structures. 

Such a large transcend regime could be identify with 
the presence of a stable QSS of the Vlasov dynamics ||7J), 
which actually evolves in the collisional timescale r cr 3> 
To. This fact has been shown throughout extensive calcu- 
lation of the average of the nonequilibrium temperature 
T D with N = 10 2 - 10 4 rotators over ~ 10 3 - 10 4 realiza- 
tions of microscopic dynamics, where it is shown agrow- 
ing of the relaxation timescale r eq with N [ll], [IH, [l3j . As 
already pointed out theoretically and observed in many 
numerical studies, the nature and the relaxation time of 
such a collisional quasi-stationary evolution depend on 
the initial conditions. 

The qualitative features of system evolution are also 
captured by the diffusion rate uj (t) — dAg (t) j dt defined 
from the square dispersion Aq (t): 

1 N 

M*) = 2AFEfoW-*<(°)] ' (10) 
which provides a measure of the average correlation func- 
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FIG. 3: Collisionless evolution of one-body distribution func- 
tion obtained from the numerical integration of the finite dif- 
ferences scheme iJHJ where it can be also appreciated the for- 
mation and destruction of filamentary structures. We illus- 
trate here the first four instants shown in the FIG[3] 
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o anomalous superdiffusion: 
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FIG. 4: Evolution of diffusion rate u (t) obtained from the 
numerical simulation of the microscopic dynamics with N = 
1000, which evidences the existence of a crossover time at 
Tcross ~ 10 4 where the anomalous superdiffusion with effec- 
tive exponent a ~ 1.68 turns a normal diffussion. 



tion: 



1 N r* 

w (*) = jz 22 (*) > with (*) = / U (*) L i ~ T ) dT > 
i=i ^° 

and accounts for the character of the diffusion regimes 
undergone by the system 1 . Dynamical studies revealed 
that after a brief ballistic regime with a = 2, the system 
exhibits an anomalous superdiffusion during the quasi- 
stationary evolution. The long time calculation of the 
diffusion rate illustrated in FIG0] allows to estimate the 
exponential constant a ~ 1.68 for the leading behavior 
of the square dispersion with N — 1000. Notice the exis- 
tence of a crossover time r cross ~ 10 4 where the diffusion 
regime turns normal, which can be related to the final 
evolution of the collisional regime towards the achieve- 
ment the of the Boltzmann-Gibbs equilibrium. Despite 
the incidence of size effects in the crossover time T cross 
and the exponent a, the general behavior of the diffu- 
sional regime described above is consistent for different 
system sizes N. 



B. Relaxation of two— body correlation function 

By denoting S [x - a* (*)] = 5 [9 - 0< (<)] S\p- Pl (t)] 
with x = (9,p) and the average over many realizations 
of the microscopic dynamics (•), we are able to introduce 



1 A regime with Aq (t) oc t a exhibits a normal diffusion when 
a = 1, ballistic diffusion a = 2, and anomalous diffusion when 
1 < a < 2 (superdiffusion) or < a < 1 (sub-diffusion). 



the one-body and the two-body distribution functions: 

/(*;*) = (^Ei=i*[*-*i (*)]), (12) 

( Xl ,X 2 ;t) = jv^jfny (E^j S [XX - Xi (t)]S[x 2 - Xj (t)] 

as well as the two-body correlation function g (xi, x 2 ; t) 
as follows: 

/ (2) (x 1 ,x 2 ;t) =f{x x ;t)f (x 2 ;t) + g (x u x 2 ; t) . (13) 

These distributions functions can be used to define the 
average (a) and the correlation function c a of a micro- 
scopic observable a (x) = a (0,p) of an individual rotator: 



(a) = I a(x)f(x;t)dx, (14) 
a(xi) a(x 2 ) g (xi,x 2 ;t) dxidx 2 , (15) 



with dx — d9dp 7 which represent two indicators of their 
corresponding dynamical evolutions. 

The correlation function of the kinetic energy ck can 
be derived from the square dispersion of the kinetic en- 
ergy per particle a 2 K = (K 2 } — (K) 2 obtained after several 
realizations of the microscopic dynamics as follows: 



ck 



N 



l —{N<T K -a 2 m } = ^^C K . 



N -I 



(16) 



where cr^ = j |(p 4 ) ~ (p 2 ) 2 } ^ s ^ ne square dispersion 
of the kinetic energy of an individual rotator, which in 
an analogous way as (K) only depends on the one-body 
distribution function. The short-time evolutions of the 
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FIG. 5: Short-time evolutions of (a) the kinetic energy (K), 
(b) the square dispersion of individual rotators and (c) 
the correlation function Ck obtained from a statistical av- 
erage over 1500 trajectories with N = 500 and 2000. The 
exponential decay of the correlation function Ck evidences 
that the two-body correlation function g (xi,X2',t) undergoes 
a relaxation process in the same timescale of the violent re- 
laxation of the one-body distribution function / (x;t). 



functions (K) , and Ck obtained from the average of 
microscopic dynamics with N — 500 and 2000 over 1500 
trajectories are displayed in FIG[5] 

While (K) and of n account for the violent relax- 
ation undergone by the one-body distribution function 
/ (9,p; t) during the first stages of the system evolution, 
the correlation function Ck evidences also the occurrence 
of a relaxation process for the two-body correlation func- 
tion g (0i, pi, 0i,Vi', t) in the same timescale, where any 
two uncorrelated particles at the initial state become cor- 
related after the incidence of such a relaxation process. 
Despite the existence of finite size effects, numerical sim- 
ulations with different N confirm the consistence of this 
observation. We shall show in a forthcoming paper that 
this result can be easily justified in terms of the well- 
known BBGKY hierarchy. Although simple, this latter 
behavior has a remarkable importance in understanding 
the collisional relaxation of this model system and the 
development of appropriate kinetic equations for the one- 
body distribution function / (0,p;t). 



C. Dynamical phase coexistence 

It is well-known that the existence of the long-range 
order leads to a very slow relaxation in the neighbor- 
hood of the critical point. However, most of works 
devoted to the numerical study of dynamical aspects 
of the HMF model also revealed the existence of very 
large relaxation times during the supcrdifussional regime 
within the energetic region u\ = 0.5 < u < u c = 0.75 
[1 H B 0, 0, IE Gj, 03 • It is remarkable that the 
above energetic region is precisely the same one where the 
microcanonical susceptibility Xm (u; b = 0) experiences a 
large increasing below of the critical point u c (see in 
FIG. 2 of the previous paper This observation sug- 

gests that the microcanonical results and these anoma- 
lous dynamical behaviors could be closely related. 

While the thermodynamical characterization of this 
model shows the features of the typical second-order 
phase transitions, from the viewpoint of its out-of- 
equilibrium dynamics, the HMF model also shows certain 
analogy with the phase coexistence phenomenon. The 
present interpretation is naturally arisen from character- 
ization of the rotator individual dynamics as a mathe- 
matical pendulum motion: when the pendulum energy 
e(0,p) — ^p 2 — mcos0 < m, the particle exhibits an 
oscillatory regime around the equilibrium configuration, 
while for e (9,p) > m the particle performs rotations in 
a given direction. The first kind of motion character- 
izes the particles trapped in the cluster existing when 
u < u c , while the second one characterizes those par- 
ticles that are not bonded to the cluster, which is the 
typical motion predominating in the homogeneous phase 
with u > u c . We can rephrase the microscopic dynam- 
ics by considering that all those rotators satisfying the 
condition e (0, p) < m belong to a dynamical clustered 
phase, while the other with e(0,p) > m belong to a dy- 
namical gaseous phase. Such a reinterpretation suggests 
the occurrence of a dynamical phase coexistence within 
the energetic region ui < u < u c . 

The relative population p of the clustered phase can 
be obtained from the distribution function f(0,p;t) as 
follows: 

/>oo 

l-p(t)= / deF(e;m,t), (17) 

J m 

where F (e; m, t) is the energy distribution function: 
F (e; m,t)= d9 d P 6 {e - e (0,p)} f (0,p; t) . 

JO J -oo 

(18) 

Since during collisional regime the distribution function 
exhibits a quasi-stationary evolution throughout stable 
QSSs of Vlasov equation ([7]), the distribution function 
can be taken as /qss (9,p;t) ~ / [e (0,p) ; m, t], which 
leads to the following form of the quasi-stationary energy 
distribution function: 



F (e; m,t)=T (e, m) f (e; m, t) . 



(19) 
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FIG. 6: Energy dependence of the relative population p of 
the clustered phase and the magnetization m (added for com- 
parative purposes). Inserted graph: the x dependence of the 
function $ (a;). This figure shows that the relative population 
of the clustered phase exhibits its most significant variation 
within the energy interval iti < u < u c . 
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FIG. 7: "Particle interchange" between the clustered and the 
gaseous phases. Panel a: Transition between vibrational and 
rotational motions during the quasi-stationary regime; Panel 
b: Evolution of the energy e (t) of a given particle and the 
magnetization m(t); Panels c and d: Frecuency counts and 
correlation function of the rate of change of the particle energy 
7 (t) = de (t) jdt. 



being T (e, m) the time period of the mathematical pen- 
dulum with energy e: 



T (e, m) 



dd 



:e,m) 



d6 



m cos ( 



(20) 

where 9 m is the positive turning point where 
p (9 m ; e, m) — when £ < m, or 9 m = ir when s > m. 

The presence of the time period function T{e 1 m) 
in Eq.(|19p demonstrates the existence of a pole in 
quasi-stationary energy distribution function at e = m 
(lim e — > m T (e, m) = oo), which separates the energetic 
range of the system from the dynamical viewpoint in a 
clustered phase when e < m and a gaseous phase for 
e > m. The incidence of such a lost of analyticity on the 
quasi-stationary regime is much significant for energies 
u < u c close to the critical point. The relatively large 
collisional relaxation times and the anomalous superdif- 
fusional observed in the interval u\ < u < u c should 
be related in some way to this dynamical phase coexis- 
tence. Unfortunately, a complete study of these question 
demands the knowledge of a kinetic equation describing 
the collisional evolution of the system. 

Nevertheless, it could be useful to calculate the relative 
population p of the clustered phase during the thermody- 
namic equilibrium of the HMF model. The integration of 
(|17p by using the Boltzmann-Gibbs distribution function 
f B G (e) oc exp (-/3e) yields: 



p(x) 



exp (x) $ (x) 



(21) 



where x — (3m, Iq (x) the modified Bessel function of 
zero-order and $ (x) is given by: 



$(x) 



exp 



y 



y- 



(22) 



being K (x) the complete Legendre elliptic integral of the 
first kind: 



K(x) 



dt 



tt/2 

yl — x 2 sin 2 1 



(23) 



The energetic dependence of the relative population 
p of the clustered phase is illustrated in FIGE1 As ex- 
pected, the system experiences an abrupt formation of 
the clustered phase below of the critical point u c . Inter- 
estingly, the transformation of the clustered phase into 
the gaseous one mainly takes place on the anomalous en- 
ergy interval u\ < u < u c , a result that reinforces the hy- 
pothesis about the origin of anomalous behavior in that 
region should be related to the occurrence of a dynamical 
phase coexistence during the out-of-equilibrium regime. 

Panel a) in FIG [7] shows that a given rotator randomly 
changes the character of its motion during the quasi- 
stationary regime, a phenomenon which could be inter- 
preted here as a " particle interchange" between the clus- 
tered and the gaseous phases. A better understanding 
of this process is achieved by analyzing the evolution of 
the particle energy e = e (9,p) = \p 2 — m • m (9) and 
the magnetization m = |m| illustrate in panel b) of this 
figure. Transitions towards the gaseous phase take place 
wherever the dependence e (t) overcome the magnetiza- 



G 



tion m (t). Although the dynamical evolution of magne- 
tization and the particle energy are actually smooth at 
the microscopic timescale To, these evolution look very 
rough at a larger timescale. It is very interesting to re- 
mark that the evolution of rate of change of the particle 
energy 7 (t) — de (t) jdt seems to be modeled as a Gaus- 
sian noise with a correlation time t c comparable to the 
microscopic timescale To, as shown in panels c) and d). 



III. PHENOMENOLOGICAL APPROACH OF 
COLLISIONAL RELAXATION 

Recently, some authors are carrying out some pioneer- 
ing studies in order to obtain some kinetic equations able 
to explain how this model system evolutes towards the fi- 
nal equilibrium configurations [lj, [H, [H, ELzi Ell ■ While 
the early relaxation dynamics is very-well described by 
the Vlasov dynamics (O, no one of such attempts sat- 
isfactorily explains the dynamical behavior of the HMF 
model in during the collisional regime. 

This aim could be achieved by specifying the dynami- 
cal mechanism allowing the energetic interchange among 
the system particles during the collisional regime. As al- 
ready commented in the previous section, the energetic 
interchange among the rotators during the collisional 
regime takes place as consequence of small oscillations or 
fluctuations of the magnetization vector m around the 
mean value (m). The characteristic timescale of such 
fluctuations is comparable to the microscopic timescale 
To, which is the characteristic timescale of the rotator in- 
dividual evolution. This latter observation and the form 
of Eq. ([3]) support the idea that the energy interchange in 
the HMF model takes place by means of the parametric 
resonance 

The mechanism of parametric resonance potentially 
appears whenever exist an oscillatory dependence of the 
microscopic parameters of an oscillatory system (2lj . 
According with the common understanding, this phe- 
nomenon can manifest in the dynamics © when a 
Fourier mode u> of the magnetization m (t) = (m) + 
a w (t) cos (tut + 5^) belongs to a very small neighbor- 
hood of the frequency Ui — 2ir/Ti of a given rotator: 
\oj - u>i\ < (a^l/wo ~ l/VN, being ui Q ~ y/(m). 

At low energies, most of the system rotators exhibit a 
harmonic oscillatory motion whose frequencies are very 
close to luq. Under these conditions the parametric reso- 
nance allows an affective energy interchange among the 
particles leading the system to a fast equilibration, which 
explains the concordance of thermodynamical predictions 
with the microcanonical numerical experiments at low 
energies @ . The growing of the energy per particle leads 
to the increasing of the non harmonic character of the 
rotator motion, and consequently, the spreading of the 
Fourier modes of m(t). Thus, the necessary conditions 
for the parametric resonance are only eventually satis- 
fied, allowing in this way to a given rotator to gain or 
lost a certain amount of energy, sometimes significant, in 



a brief time period, as already illustrated in FIG|7l Such 
large energetic transfers are rather analogous to a close 
encounter. 

Recently studies about the origin of the Hamiltonian 
chaos by using the methods derived from the Rieman- 
nian interpretation of dynamics conclude that the phe- 
nomenon of parametric resonance is the fundamental 
mechanism to induce irregular motion in most of the 
Hamiltonian systems which are binding for most of re- 
gions of the configurational space [12, HE HH HH ■ Ap- 
parently, the parametric resonance is the dynamics me- 
chanics explaining the origin of dynamics instability and 
the collisional relaxation of the HMF model. 

When the dynamical behavior of the system param- 
eters is very complex, the treatment of such a micro- 
scopic dynamics could be performed by using the meth- 
ods of theory of random flights (see in ref. [24j | the appli- 
cation of such methods in the framework of the Hamil- 
tonian chaos). This idea follows from the observation 
that the small oscillatory behavior of the magnetiza- 
tion vector in timescale longer than the microscopic 
time To may be taken into account as a Gaussian noise. 
In fact, the dynamical evolution of the rotator energy 
e = \p 2 — m • m (6) illustrated in FIG [7] suggests that 
the rotator dynamics is rather analogous to a Brownian 
motion, that is, a diffusion-like process. 

Let us apply the present hypothesis in order to obtain 
a phcnomenological kinetic equation for describing the 
collisional evolution of the HMF model. Hereafter, the 
Einstein summation convention is assumed. The magne- 
tization vector m could be decomposed by considering 
the its mean value (m) plus a Gaussian noise term as 
follows: 

m ~ (m) + n a ry a , (24) 

where the unitary vectors n a consider a different behav- 
ior of the noises n a in the parallel and transverse direc- 
tion of the mean value (m) . Let us also assume that the 
noise only depends on the global system observables (sys- 
tem size N, energy u, average magnetization vector (m), 
etc.). Since a purely diffusive influence does not lead 
the one-body distribution function towards the Gaussian 
profiles of velocities resulting at the final equilibrium, we 
shall take into consideration the inclusion of a dynamical 
friction in this context in order to accomplish this feature 

Thus, we conjecture that the microscopic dynamics of a 
given rotator could be conveniently approximated by the 
following Langevin equation with a multiplicative noise: 

9 = p, p = -\(9)p + m(9) x (m) + e a (9)n a . (25) 

where the dynamical friction A (9) should necessarily be 
^-dependent (see in Eq. (j30|) below). The functions e a (8) 
are derived from the unitary vectors n a as e a (9) = 
m (8) x n a . The Gaussian noises n a satisfy the following 
relations: 
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(%(*)> = 0, {r] a (t)r) b [t')) = 2n ab 6{t-t'). (26) 

where r2 t is the correlation matrix. The Langevin equa- 
tion (|2"5|) leads to the following Fokker-Planck equation: 



df df df d 



df 
dp 



VPf 



(27) 

where the mean value of the magnetization vector (m) 
was taken as follows: 



<m)=m[/] = J d9dpm[9) /, 
and the diffusion coefficient D [9] by: 



D [9] = n ab e a (9) e b (9) 



(28) 



(29) 



We also assume that D [9] is related to the dynamical 
friction by means of the Einstein relation: 



A (6>) = j]D (9) 



(30) 



being rj the nonequilibrium analogue of the inverse tem- 
perature, rf 1 = J p 2 f (9,p;t) dpd9. The consistence of 
Fokker-Planck equation (|27[) demands that the correla- 
tion matrix £l ab should be considered as certain func- 
tionals of the distribution function /. Since il ab de- 
creases with the increasing of the system size, the dif- 
fusive term of Fokker-Planck equation is only effective at 
large temporal scales. Thus, the assumption (|28p takes 
into account that the dynamical evolution of the distri- 
bution function at the timescale To is described by the 
collisionless dynamics (|7|), while diffusive term considers 
the relaxation towards the Boltzmann-Gibbs distribution 
function at the collisional timescale r cr oc D~ l . 

The definition of the diffusion coefficient D (9) can be 
conveniently rewritten by considering the tridimensional 
vector m = (cosfl, sinfl, 0). Thus, the scalar term e a e b 
could be rephrased as: 



e a e b = (m x n ) ■ (m x n b ) = 5 ab - (m • n a ) (m • n Q ) . 

(31) 

Let us now reconsider again the bidimcnsional form of 
the vector m (9) = (cos 9, sin 9). By denoting the z-th 
component of m (9) by (9) , the diffusion coefficient 
can be expressed as: 



D = m,i (9) K-ijUij {9) , 



(32) 



where ICi 



Q5t 



tt ab {n a )i [nb) a, and fl = Sp [Q ab 



mean value of the magnetization (m), the symmetric ma- 
trix ICij should take the form: 

fcij = a [/] % + b [/] [mi) (mj) , (33) 
allowing to rewrite the diffusion factor as follows: 



D(9)=a[.f}+b[f}{ m [f}-m(9)Y 



(34) 



Since the only preferential direction of this system is the 



The closure of the Fokker-Planck approximation relies on 
the determination of the functional a [/] and b [f] , which 
is a task beyond of our phcnomcnological reasonings. It 
is very easy to verify that Eq. ([2"T]) ensures the particles, 
momentum and the energy conservation. 

Chavanis et al introduced in r efs.0, EE [3 a sim- 
ilar treatment by modifying of the HMF model micro- 
scopic dynamics ^ with the inclusion of frictional and 
stochastic terms, whose resulting Fokker-Planck equa- 
tion basically differs from the present proposal by the 9- 
dependence of the diffusion coefficient D (9). The origin 
of this difference is found in the fact that these authors 
actually introduce a variant of HMF model (Brownian 
Mean Field model) which accounts for the external influ- 
ence of a thermostat, while the frictional and stochastic 
terms considered in the present work try to capture the 
effective influence of the magnetization fluctuations on 
the exact individual dynamics of a given rotator in a 
timescale larger than the microscopic time To. The cru- 
cial question now is how to justify and precise within the 
framework of the parametric resonance the ad hoc ap- 
proximation (|25[) . mainly, the inclusion of the effective 
dynamical friction. 



IV. FINAL REMARKS 

As already appreciated in many studies, the HMF 
model dynamics undergoes the incidence of two relax- 
ation process with different characteristic timescales: (1) 
the violent relaxation timescale T vr leading to the quasi- 
stationary evolution with a superdiffusion regime, and 
(2) the collisional relaxation timescale r cr where the 
Boltzmann-Gibbs equilibrium is achieved as final state. 
It is well-established theoretically and numerically that 
the collisionless relaxation timescale is the same order of 
the microscopic timescale, r vr ~ To, while the collisional 
relaxation timescale seems to grow with the system size 
in some power of N, r cr ~ tqN". Early estimations pre- 
dict a simple linear dependence 0], although there are 
also several numerical ex per iments evidencing a nontriv- 
ial exponent a ~ 1.7 [II []J, [l3[ . 

Generally speaking, the determination of the correct 
A^-dependence of the collisional relaxation time is rather 
difficult to carry out by means of numerical compu- 
tations of the microscopic dynamics since such exper- 
iments are plagued of prejudicial size effects and poor 
equilibration of trajectory averages whose elimination 
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seems to demand the consideration of more larger num- 
ber of particles than the actual computational possibil- 
ities N ~ 10 4 — 10 5 For example, besides of the 
nontrivial exponent a = 1.7 of the power-law growing of 
r cr reported by Yamaguchi et al in refs.[L^[l3|. they also 
indicated the following A^-dependence of the quantity q: 



d(m (t)) 



d(\nt) 



oc VN 



(35) 



If the above estimation of r cr is the correct relaxation 
timescale associated to the collisional quasi-stationary 
evolution (qse) of the one-body distribution function 
/ (6,p;t), the temporal dependence for large N should 
be given by: 



/ (0, p; t) ~ f qse (9, p; t/r cr ) +0(1/N). 



(36) 



The above assumption implies that the leading behavior 
the temporal dependence of the average magnetization 
should be given by (m(t)) ~ A(t/r cr ), where A (s) — 
|J m (9) fo (0,p; s) d6dp\ is a size-independent function. 
Thus, the leading behavior of the quantity: 



q = 



d(m (t)) 



d\nt 



_ dA(s = 1) 
ds 



(37) 



should be N-independ when N is large enough. Since the 
numerical result (|35[) disagrees with this reasoning, the 
nontrivial exponent a ~ 1.7 reported in the Yamaguchi 
et al study might not correspond to the true scaling be- 
havior for N very large, but only a finite size behavior 
which should disappear for N large enough. 

The size dependence of the collisional relaxation 
timescale is a question with a primordial importance in 
order to understand the dynamical features of a long- 
range interacting system as the HMF model when ther- 
modynamic limit is invoked. As already pointed out 
by some authors, the existence of relaxation timescale 
T cr diverging with the imposition of the thermodynamic 
limit, lirn/v^oo T cr = oo, leads automatically to the 
non commutative character of the the thermodynamic 
limit limjv^oo with the infinite time limit lim^oo neces- 
sary to the equilibration of temporal averages (Azv) t = 
f* A N (r) dr/t: 

ton lim iMi ± um i im <Mi, (38) 



Thus, the imposition of the thermodynamic limit before 
the infinite time limit makes endless the duration of col- 
lisionless dynamics, and therefore, the Boltzmann-Gibbs 



equilibrium arises as an admissible but particular quasi- 
stationary state of the violent relaxation. Such a dy- 
namical picture is not only a feature of the HMF model 
dynamics. On the contrary, the large the collisional re- 
laxation times estimates of many astrophysical objects 
like the elliptical galaxies (r cr ~ OAtqN/ In A" where mi- 
croscopic timescale to ~ 1/ y/TTp) induces also to suppose 
the collisionless character of the dynamical evolution of 
such real long-range interacting systems, which explains 
the rich variety of structures observed in this context [2(| • 

Let us use the phenomenological approach of collisional 
regime in terms of the Fokker-Planck equation ([77)) in or- 
der to analyze the Af-dependence of the collisional relax- 
ation timescale r cr . Since the amplitude of the Gaussian 
noises rj a describing the magnetization fluctuations dur- 
ing the quasi-stationary evolution decreases as 1/ y/N as 
well as the underlying correlation times are of order of 
microscopic timescale To (see in panel c) of FIG|7]), the 
correlation matrix f2 a f, (|26[) and the diffusion coefficient 
D (9) ([25]) decrease as l/N. Thus, the phenomenological 
picture described in this work suggests that the relevant 
timescale of the collisional evolution of the HMF model 
should obey a linear growing r cr = t^N . 

As already discussed in our previous work [li|, such a 
linear A^-dependence of r cr seems to be intimately related 
to the Af-dependence of the additive constant of entropy 
per particle sq = \ In (2ire 2 Ig/N^ , since the imposition 
of a scaling dependence of the coupling constant g as 
g (N) — instead of the usual Kac prescription [27j 
gkac (N) = 7 /N leads to a simultaneous regularization 
of the divergence of the collisional relaxation timescale 
T cr = tqN = \fTWJg = y/Tpy and the additive constant 
,s = \ In (2we 2 Ig/N) = \ In (27re 2 /7) in the thermody- 
namic limit , avoiding in this way the incidence of an 
undesirable dynamical anomalies like the one described 
in Eg. (|3"8]) and the divergence of the thermodynamic po- 
tentials. 

Although the phenomenological approach of collisional 
relaxation described in this work supports the validity of 
the linear growing r cr = tqN, the problem about the 
exact A^-dependence of collisional relaxation timescale 
could be definitively solved by means of the development 
of appropriate kinetic equations starting from the con- 
sideration of first principles [l6| . We shall address in our 
forthcoming paper an intense investigation of the HMF 
model dynamics in terms of kinetic equations, where 
the implementation of suitable kinetic equations start- 
ing from the well-known BBGKY hierarchy will receive 
a primordial attention. 
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